Numerical Simulations of Vortex Shedding in Hydraulic Turbines 


Daniel Domey* 

NASA George C. Marshall Space Flight Center 
Marshall Space Flight Center, AL 

Bogdan Marcu* 

The Boeing Company 
Rocketdyne Propulsion and Power 
Canoga Park, CA 


ABSTRACT 

Turbomachines for rocket propulsion applications operate 
with many different working fluids and flow conditions. 
Oxidizer boost turbines often operate in liquid oxygen, 
resulting in an incompressible flow field. Vortex shedding 
from airfoils in this flow environment can have adverse 
effects on both turbine performance and durability. In this 
study the effects of vortex shedding in a low-pressure 
oxidizer turbine are investigated. Benchmark results are 
also presented for vortex shedding behind a circular 
cylinder. The predicted results are compared with 
available experimental data. 

INTRODUCTION 

Modern high-work turbines can be compact, transonic, 
supersonic, counter rotating, and can use a dense drive 
gas. The vast majority of modem rocket turbine 
designs fall into these categories. These turbines are 
often characterized by large amounts of flow 
unsteadiness. The flow unsteadiness can have a major 
impact on turbine performance and durability. For 
example, the Space Transportation Main Engine 
(STME) fuel turbine, a high- work transonic design, 
was found to have an unsteady interrow shock which 
reduced efficiency by 2 points and increased dynamic 
loading by 24 percent. The Revolutionary Reusable 
Technology Turbopump (RRTT), which uses full flow 
oxygen for its drive gas, was found to shed vortices 
with such energy as to raise serious blade durability 
concerns. In both cases, the sources of the problems 
were uncovered before turbopump testing with the 
application of validated unsteady computational fluid 
dynamics (CFD) to the designs. In the case of the 
RRTT and the Alternate Turbopump Development 
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(ATD) turbines, the unsteady CFD codes were used not 
just to identify problems, but also to guide designs that 
mitigate problems due to unsteadiness. Using unsteady 
flow analyses as part of the design process has led to 
turbine designs with higher performance and fewer 
dynamics problems. References 1-4 are examples of 
the application of unsteady CFD to rocket turbine 
designs. 

The effects of vortex shedding in a low-pressure 
oxidizer turbine are focus of the current study. The 
numerical analysis is first tested on the benchmark case 
of a cylinder in water, and is then applied to the turbine 
operating in liquid oxygen. The predicted results have 
been compared with the available experimental data. 

NUMERICAL ANALYSIS 

The governing equations in the code, called 
AARDVARK, are the quasi-three-dimensional, 
unsteady, Navier-Stokes equations. The equations have 
been written in the Generalized Equation Set (GES) 
format (a detailed description of the GES can be found 
in Ref. 5). A description of the AARDVARK 
code/algorithm, as well as its application to several 
turbine and pump test cases, is presented in Refs. 6 and 
7. The general code structure is based on a well- 
established compressible, quasi-three-dimensional, 
unsteady turbomachinery flow solver [8]. It employs a 
system of overset O-grids and H-grids, with the values 
on the O-H boundaries being updated each time step by 
bilinear interpolation from the adjacent grid. The 
Baldwin-Lomax turbulence model is used for 
turbulence closure [9]. The code contains two options 
for the fluid properties. The first option is based on the 
equations of state, thermodynamic departure functions 
and corresponding states principle constructed by 
Oefelein [10]. The second option, used only for liquids, 
is based on splines generated from the NIST Tables 
[ 11 ]. 
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CYLINDER BENCHMARK SIMULATION 

The grid topology for the cylinder contains 361x71 
points in the O-grid and 142x121 points in the H-grid, 
for a total of 42,813 points (see Fig. 1). The average 
value of y+, the non-dimensional distance of the first 
grid point above the surface, is approximately 0.1. 



Figure 1. Computational grid for the cylinder. 

The Reynolds number for the benchmark simulation, 
Re=30,000, was chosen to approximate the Reynolds 
number based on the oxidizer turbine first-stage vane 
trailing edge diameter. 

The vortex street behind the cylinder is highlighted 
using non-dimensional velocity contours in Fig. 2. The 
vortex street is observed to convect through the exit 
boundary with little or no reflection. 



Figure 2. Non-dimensional velocity contours about 
the cylinder. 


The vortex shedding induces a periodic variation in the 
static pressure along the back half of the cylinder. 
Figure 3 shows the variation of the non-dimensional 
static pressure as a function of time near the trailing 
edge stagnation point on the cylinder. 



Figure 3. Unsteady static pressure near the trailing 
edge stagnation point. 

Performing a Fourier decomposition on the unsteady 
static pressure yields a Strouhal number of St=0.23. 
The predicted Strouhal number is greater than the 
experimentally observed value of St=0.20 [12], but still 
considered accurate within the limitations of the 
turbulence modeling and grid resolution. 

LOW-PRESSURE OXIDIZER TURBINE STAGE 

The turbine stage under consideration is the first-stage 
of the Space Shuttle Main Engine (SSME) Low 
Pressure Oxidizer Turbine (LPOT). The SSME LPOT 
was recently the subject of an in-depth investigation 
because cracks were discovered on the first stage vanes 
[13]. Computational simulations were used to 
determine that vortex shedding was the probable cause 
of the cracks, and experimental data was used to 
pinpoint the shedding frequency at approximately 36 
kHz. 

The grid topology for each vane contains 301x71 points 
in the O-grid and 152x121 points in the H-grid. Each 
rotor contains 301x71 points in the O-grid and 135x81 
points in the H-grids. The 2-vane/3-rotor simulation 
contains a total of 176,444 grid points (see Fig. 4). The 
average value of y+ is approximately 0.5 for both the 
vane and the rotor, and the boundary layers are 
discretized with an average of 25 grid points. 

The operating fluid is liquid oxygen (LOX) and the 
pressure ratio across the stage is approximately 0.79. 
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Figure 4. Computational grid for the SSME LPOT. 

Instantaneous and time-averaged non-dimensional static 
pressure contours are shown in Figs. 5 and 6, 
respectively, while instantaneous and time-averaged 
non-dimensional velocity contours are shown in Figs. 7 
and 8. 



Figure 5. Instantaneous non-dimensional static 
pressure in the SSME LPOT turbine. 



Figure 6. Time-averaged non-dimensional static 
pressure in the SSME LPOT. 


The instantaneous static pressure contours hint at vortex 
shedding emanating from the trailing edges of both the 
vane and rotor, but do not identify it in the interblade 
region. 



Figure 7. Instantaneous non-dimensional velocity in 
the SSME LPOT. 
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Figure 8. Time-averaged non-dimensional velocity in 
the SSME LPOT. 


As expected the time-averaged pressure exhibits no 
traces of vortex shedding (see Fig. 6). The 
instantaneous velocity contours indicate significant 
vortex shedding from both the vane and rotor airfoils. 
The vortex shedding from the vane remains coherent all 
the way into the rotor passage. Similar to the time- 
averaged static pressure, the time-averaged velocity 
contours do not reveal the presence of the vortex 
shedding. 

Time-averaged pressure distributions along the airfoil 
surfaces are shown in Fig. 9. 



x/c 


Figure 9. Time-averaged surface pressure 
distributions. 

The vane airfoils are aft-loaded, while the rotor blades 
are relatively evenly loaded. In addition, the majority 
of the pressure drop occurs across the vane. Figures 10 
and 1 1 contain an unsteady pressure trace at the vane 


trailing edge and the corresponding Fourier 
decomposition. 



Figure 10. Unsteady pressure trace at the vane 
trailing edge. 



Figure 11. Fourier decomposition of the unsteady 
pressure trace at the vane trailing edge. 

The uncorrected predicted frequency for the vortex 
shedding is higher than that deduced experimentally 
(42.5 kHz vs 36 kHz). Note, however, that the 
benchmark simulation for the circular cylinder at the 
same Reynolds number yielded a Strouhal number that 
was approximately 15% too high. If this benchmark 
correction is applied to the SSME LPOT, then a 
frequency of 36.9 kHz is obtained. 

The final paper will contain a more in-depth 
discussion of the results. 
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CONCLUSIONS 

Unsteady numerical simulations have been performed 
for a circular cylinder operating in water and a low- 
pressure oxidizer turbine operating in liquid oxygen. 
The circular cylinder test case was used to 
validate/anchor the analysis for unsteady flows in 
liquids. The oxidizer turbine test case was used to study 
vortex shedding in a turbine that had developed cracks 
on the first-stage vanes. The predicted shedding 
frequency for the vanes, after applying a correction 
factor based on the cylinder simulation, is in agreement 
with the available experimental data 

Additional conclusions will be added based on further 
interrogation of the results . 
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